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^ , A brief review of recent developments in the theory of the Riemann zeta function inspired 

• by ideas and methods of quantum chaos is given. 

o ; 

en , §1. Introduction 

Q ■ At the first glance the Riemann zeta function and quantum chaos are completely 

, disjoint fields. The Riemann zeta function is a part of pure number theory but 

quantum chaos is a branch of theoretical physics devoted to the investigation of 
non-integrable quantum problems like the hydrogen atom in external fields. 
^ , Nevertheless for a long time it was understood that there exist multiple inter- 

I relations between these two subjects. ^^'^^ In Sections [2] and [3] the Riemann and the 

I ■ Selberg zeta functions and their trace formulae are informally compared. '^^ From 

, the comparison it appears that in many aspects zeros of the Riemann zeta function 

I resemble eigenvalues of an unknown quantum chaotic Hamiltonian. 

CN ■ One of the principal tools in quantum chaos is the investigation of statistical 

, properties of deterministic energy levels of a given Hamiltonian. In such approach 

one stresses not precise values of physical quantities but their statistics by con- 
' sidering them as different realizations of a statistical ensemble. According to the 

! BGS conjecture"^^ energy levels of chaotic quantum systems have the same statistical 

' properties as eigenvalues of standard random matrix ensembles depended only on 

the exact symmetries. In Section U] it is argued that is quite natural to conjecture 
that statistical properties of the Riemann zeros are the same as of eigenvalues of the 
' gaussian unitary ensemble of random matrices (GUE). This conjecture is very well 

confirmed by numerics but only partial rigorous results are available.^-* 

In Section [5] a semiclassical method which permits, in principle, to calculate 
correlation functions is shortly discussed. The main problem here is to control cor- 
relations between periodic orbits with almost the same lengths. In Sections [6] and 
[7] it is demonstrated how the Hardy-Littlewood conjecture^^ about distribution of 
near-by primes leads to explicit formula for the two-point correlation function of 
the Riemann zeros. The resulting formula describes non-universal approach to the 
GUE result in excellent agreement with numerical results. 

In Section [8] it is demonstrated how to calculate non-universal corrections to the 
nearest-neighbor distribution for the Riemann zeros. 

Spectral statistics is not the only interesting statistical characteristics of zeta 
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functions. The mean moments of the Riemann zeta function along the critical line 
is another important subject that attracts wide attention in number theory during 
a long time. In Section [9] it is explained how random matrix theory permit Keating 
and Snaith to propose the breakthrough conjecture about mean moments. This 
conjecture now is widely accepted and is generalized for different zeta and L-functions 
and different quantities as well.^''-*'^^^ 



§2. Riemann and Selberg zeta functions 



The Riemann zeta function, Ci^), is defined as 

oo ^ 

C(-) = E;^ (2-1) 

n=l 

where the summation is taken over all integers and s is a complex parameter (see 
e.g.^^-'). As each integer can uniquely (up to the ordering) be written as a product 
over prime numbers p, n = p^^p^^ ■ ■ ■ p^*" , the sum in (|2-ip may be expressed as the 
Euler product over all prime numbers 

The Selberg zeta function, Z{s), is related with hyperbolic motion on constant cur- 
vature surfaces generated by discrete groups. Similarly to (j2-2p it is defined as 
the product but not over primes but over all primitive periodic orbits (ppo) for the 
motion on the surface considered 

oo 

Z(s) = Jl ]J(l-e-^^'(^+™)) . (2-3) 

ppo m=0 

Here Ip are the lengths of these orbits and s is a complex number. (An informal 
introduction to this subject is given e.g. in Ref. 14).) 

The zeta functions (j2-2p and (|2-3p converge only when Re(s) > 1 but the both 
can analytically be continued into the whole complex plane of s where they obey 
similar functional equations 

C(s) = 0r(s)C(1-s), Z{s) = <Ps{s)Z{1-s) (2-4) 

where for the Riemann zeta function 

0i?(s) = 2*^*-^ sin (7rs/2) r(l - s) . (2-5) 

For the Selberg case the function (/>s(s) depends on the group considered. In the 
simplest case of groups with compact fundamendal domain with area fj, 

fS-l/2 \ 

4>s{s) = exp I fi I titanvrn du . (2-6) 
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The both zeta functions have also a striking resemblance of their zero structure (see 
Fig. [T|). There exist two types of zeros: trivial ones coming from known zeros of 
the factor (j)(s) and non-trivial ones which lie on the line Re(s) = 1/2. For the 
Selberg zeros this property is a theorem because its non-trivial zeros are related 
with eigenvalues of the self-adjoint Laplace-Beltrami operator on constant negative 
curvature surfaces. For the Riemann zeros it is the content of the famous Riemann 
conjecture which is confirmed by all existing numerical calculations. 



trivial zeros 



-4 -3 -2 -1 



Re s=l/2 



>/ non-trivial 



trivial zeros 



\ 



-4 -3 -2 -1 1 



Re s=l/2 



.. non-trivial 



Fig. 1. Structure of zeros for the Riemann (left) and the Selberg (right) zeta functions. The Rie- 
mann zeta function has the pole at s = 1 and trivial zeros are at s = —2, —4, . . . ,. Trivial zeros 
for the Selberg zeta function with compact domain are at s = 1, 0, —1, —2, . . . ,. 



§3. Trace formulae 

Non-trivial zeros are not known analytically (cf. Ref. 15) where lowest Riemann 
zeros are given with 1000 digits) but for the both zeta functions there exist trace 
formulae which express a sum over all non-trivial zeros, Sfi^ through a sum over 
primes for the Riemann case and over periodic orbits for the Selberg case (see e.g. 
Ref. 3)). 

Let h{r) be a test function with the following properties 

• h{r) is analytical in |Im r| < 1/2 + 5, 

• h{—r) = h(r), 

• \h{r)\ <A{l + \r\)-^-\ 

and g{u) be the Fourier transform of h{r) 

The trace formula for Riemann zeros (Weil explicit formula) is 

E h{sn) = ^J ^ Hr)jT{-^ + ^r)dr - (3-1) 

non-trivial 



zeros 
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oo ^ . . 

-2 J2 lnpY.-;^g{nlnp)+hC^) + h{-'-)-gmn7r 

primes n=l 

The trace formula for zeros of the Selberg zeta function (for groups with finite fun- 
damental domains) has similar form 

E Ksn) = ^ h{r)r tanh(vrr)dr + E E 2sinh(n/,/2) ^("^^^ ' 

(3-2) 

Here Ip are length of periodic orbits and ^ is the hyperbolic area of the fundamental 
domain. 

The Selberg trace formula is exact. For general chaotic systems Gutzwiller^^^'^'''' 
obtained semiclassical trace formula which, in general, is only the first dominant 
term in the semiclassical limit /i ^ (though the calculation of correction terms 
are in principle possible^^-*). For a two dimensional chaotic billiard with area // the 
Gutzwiller trace formula for the density of states in 'semiclassical' limit — > oo is 
the sum of two terms 

d{E) = Y,^{E- En) = d{E) + (E) (3-3) 

n 

where the smooth part is given by the usual Thomas-Fermi expression: d{E) = /u/47r, 
and the oscillatory contribution is represented as the sum over all classical periodic 
orbits with coefficients obtained from purely classical mechanics 

rp CO 1 

^{osc) (^E) = y^y cos ( n{klp - -/x„) ) . (3-4) 

ppo n=l ' ^ P 'I 



Here A; = yi? is the momentum and for each primitive periodic orbit, Ip is its 
length, Mp is the monodromy matrix, and /ip is the Maslov index. "'^'^^ 

Ignoring questions of convergence, trace formulae for the Riemann and the Sel- 
berg zeta functions, (|3-ip and (j3-2p . can also be rewritten in similar form by taking 
/i(r) = (5(r — E) and conserving only terms dominant at large energy. In particular 
for the Riemann zeros such 'physical' trace formula takes the form (13-30 with 



and 



J(E) = ±1„| (3.5) 



1 CO , ^ CO ^ 



TT — ' ^ — • p"/^ vr ^ — • ^/n 

p n=l '■ n=l ^ 



where A{n) is von Mangoldt function 



^(^) \ otherwise ' ^'^ 



Riemann zeta function and quantum chaos 



5 



The above formulae demonstrate a striking resemblance between the Riemann and 
the Selberg zeta functions. Different quantities for number-theoretical zeta functions 
have analogs in dynamical zeta functions (and vise verso) according to the following 
dictionary^)' ^) 



Riemann zeta function 


Selberg zeta function 


Primes — > 


< — Periodic orbits 


Inp — > 


i — Period Tp 


non-trivial zeros — > 


< — eigen momenta 


Prime number theorem: 


Periodic orbit number theorem: 







It opens an exiting possibility (cf. Refs. 1),3)) that the Riemann zeta function 
is similar to a dynamical zeta function for a chaotic quantum system corresponding 
to the motion on a constant negative curvature surface. Till now all attempts to 
find such system failed (cf. Refs. 19) and 20)) but a quete of a self-adjoint operator 
whose eigenvalues are related with non-trivial Riemann zeros (as it seems to be 
suggested by Polya and Hilbert) continue. Nevertheless notice that the overall signs 
of the trace formulas for the Riemann zeta and dynamical systems are opposite 
which complicates such interpretation (cf. in this connection Connes' absorbtion 
spectrum^^) ) . 

§4. Spectral statistics 

If one accepts that the Riemann zeta function is analogous to a dynamical zeta 
function and its zeros are in a certain sense related to eigenvalues of an unknown 
Hamiltonian it is tempting to applied to Riemann zeros numerous results derived 
within dynamical quantum chaos. 

First of all, it is natural to investigate statistical properties of Riemann zeros. Of 
course, Riemann zeros are defined by a deterministic procedure and are not random. 
Nevertheless, if one takes groups of N consecutive zeros in different part of the 
spectra such collection as a whole can be considered as originated from realizations 
of a random ensemble. 

According to the Bohigas, Giannoni, Schmit conjecture^) energy levels of chaotic 
quantum systems on the scale of the mean level density are distributed as eigenvalues 
of standard random matrix ensembles depended only on symmetry properties. There 
are three main classes of universality.^^) Integer spin systems without (resp. with) 
time-reversal invariance are described by Gaussian Unitary Ensemble (GUE) (resp. 
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Gaussian Orthogonal Ensemble (GOE)) and systems with half-integer spin and time- 
reversal invariance are supposed to belong to Gaussian Symplectic Ensemble (GSE). 

Symmetry properties manifest primarily in the degeneracies of periodic orbit 
lengths. Time-reversal invariant systems should have doublets of periodic orbits with 
exactly the same length, but systems without time-reversal invariance, in general, 
have no exact degeneracies. In the above analogy prime numbers play the role of 
periodic orbits. As primes are not degenerated the conjectural dynamical system 
connected with the Riemann zeros should belong to the universality class of time- 
non-invariant systems and, consequently, it should have spectral statistics as gaussian 
unitary ensemble of random matrices (GUE). 

It is well known^^) that the joint distribution of eigenvalues, Ej, for this ensemble 
has the following form 

N 

P{Ei,E2,...,EN)^l[\Ei-E,fexp{-Y,El) . (4-1) 

i<j k=l 

The knowledge of the joint distribution of eigenvalues permits formally to calculate 
n-point correlation functions which are defined as the probability density of having 
n eigenvalues at given positions 

Rn{Ei,...,En)^ j PiEi,...,EN)dEn+i...dEN . (4-2) 

For classical random matrix ensembles correlation functions are know analytically.^^) 
For GUE (j4Tp they have the determinantal form 

Rn{Ei,...,En) = det[K{E„Ej)]\ij=i,,„^n (4-3) 

where in the universal limit of unit mean density and — > oo the kernel K{Ei, Ej) 
is 

K{E,,E,) = '^^^^^^^^1^. (4-4) 
In particular, the unfolded two-point correlation function has especially simple form 

^,(e, = l-(!i^)^l-^ + ^(e-+e-) (4.5) 

where e = Ei — E2. 

The above discussion leads to the conjecture that correlation functions for high 
Riemann zeros are given by Eqs. ()4-3p - ()4-5p . It was Montgomery^) who in 1973 first 
proposed this conjecture. He considered two-point correlation function of Riemann 
zeros and proved the result equivalent to the existence of l/(27r^e^) term in ()4-5p . 
Physically Montgomery's theorem corresponds to the diagonal approximation (see 
below) . Though oscillating terms in (j4-5p cannot be obtained rigorously, Montgomery 
conjectured that spectral statistics of Riemann zeros coincides with GUE statistics. 

Odlyzko has performed very large scale computation of Riemann zeros. '^^^ In 
Fig. [2] left the two-point correlation function numerically computed by Odlyzko is 
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Fig. 2. Left: Two-point function for 2 ■ 10* Riemann zeros near the 10^^-th zero. Solid line: GUE 
prediction (|4-5|) . Right: Nearest-neighbor distribution for 10^ Riemann zeros close to i5 = 
2.5 • 10^^ Solid line: GUE prediction. 

presented. It is clearly seen that the agreement between numerical calculations and 
GUE prediction (|4-5p is excellent. Another important statistics of Riemann zeros, 
namely, the distribution of the nearest neighbor is shown in Fig. [2] right. Once more 
the agreement with GUE is striking. 

§5. Correlation functions 

In view of such excellent confirmation of GUE conjecture for Riemann zeros it 
is natural to try to derive (at least physically) correlation functions from the first 
principles. 

Formally the n-point correlation function is defined as the probability of having 
n energy levels at given positions. 

Rniei, ...,en) = {d{E + ei)d{E + es) . . . d{E + e„)) . (5-1) 

Here (...) denotes an energy smoothing 

{f{E)) = j f{E')a{E-E')dE' . (5-2) 

cf{E) is a smoothing function picked at zero (like the Gaussian). It is assumed to be 
normalized by the condition j a{E)dE = 1 and to have a certain small width AE. 

The simplest way of the using these formulae is to substitute instead of the true 
density d(E) its semiclassical expression 

diE) = d{E) + ^ Ap,„e^"^^(^)/'^ + cc. . (5-3) 

It gives the following explicit formula for the two-point correlation function 

-R2(ei,e2) = + X] ^Pi,ni^p2,n2 (^^P ^ [niSp^{E + ei) - n2Sp^{E + €2)]) + cc. 

(5-4) 
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where Sp{E + e) = Sp + TpC and T{E) is the classical period of motion. Finally one 
gets 

-R2(ei,e2) = + ^ ^Pi,ni^p2,n2-^(^'i'"i)exp ^ [niTp^{E)ei - n2Tp^{E)e2] + c.c. 

Pi,ni 

(5-5) 

and 

B{pi,ni) = (^■yi-p^[niSp^{E) - n2Sp^{E)]^ . (5-6) 

The calculation of this quantity is the main problem in the semiclassical approach 
to spectral statistics. To get the answer correctly one has to take into account tiny 
correlations between lengths of different periodic orbits which is in general a quite 
non-trivial task.^^^ 

§6. Hardy-Littlewood conjecture 



Using the trace formula for the Riemann zeta function (j3-6p one gets that the 
connected two-point correlation function of Riemann zeros -/?2 ^ — ^2 — ^ 

i?W(ei,e2) = ^ ^^^^^^ ^,i(i;+.01nn.-i(E+e.)lnn.^ ^^^^^ (g.^) 

''^ ni,n2 V^l^S 



Here (...) denotes an averaging over E as in (l5-2p . The simplest approximation 
(called the diagonal approximation^^)) consists in taking into account only pairs of 
terms with exactly the same periodic orbit length (i.e. ni = n2). Then 

Rf"'\ei,e2) = ^Y1 ^^e^(^i-^^)i°" + c.c. . (6-2) 

n 

By straightforward calculations this expression may be transformed as follows (e = 
ei - 62) 

Rf^'\e) = In (|C(1 + ie^'P^'-'He)) (6-3) 

and 

is given by the convergent sum over primes 

°° 1 _ 

^(rfia<;)(g) = exp(V y —^ei^^^P^ + c.c.) . (6-4) 

YY^^ qryf^ 

p m=l 

In the limit e ^ 0, ({I + ie) (ie)~^ and ^('^*"5)(e) const. Hence 



1 



2 y-^- 27r2e2' 

which agrees with the smooth part of the GUE result (j4-5p and corresponds to 
Montgomery's theorem^) . 
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The calculation of off-diagonal contribution is more difficult. Formally it can be 
expressed as the double sum over primes 

The exponent exp(i£^ ln(ni/n2)) oscillates quickly except for ni close to n2- One 

writes ui = n-2 + r and ln(ni/n2) ~ r/n2. In this manner R^2^^\^i-> £2) is rewritten 
as follows 

4°^^) (6) = ±,y ^^^^^^^ + "^ /ei^Vn+i.lnn\ + _ (g.g) 

47r^ n \ / 

n,r 

To proceed further one needs to control the behavior of the the product of two von 
Mangolt functions A{n)A{n + r) which is nonzero only when both n and n + r are 
power of prime numbers. 

It is naturally to assume that the dominant contribution to the above sum comes 
from the mean value of this product 

1 ^ 

air) = lim — A{n)Ain + r) . (6-7) 

n=l 

No exact theorems fix the value of this limit. It is known only from the famous Hardy- 
Littlewood conjecture^) according to which it is expresses through the following 
singular series 

Here the sum is taken over all integers q= 1,2,... and all integer p < q co-prime to 
q. ii{q) is the Mobius function 

' 1 if g = 1 
^i{q)=l {-If ifq=pi...pk (6-9) 
if g is divisible on 

and (p{q) is the Euler function which counts integers less than q and coprime to q 

p\q 

Using properties of such series one proves^) that a(r) for even r can be represented 
as the finite product 

"W=^2nF^- (6-11) 

p\r 

The product is taken over all prime divisors of r except 2 and C2 is called the twin 
prime constant 

^2 = 2 n (1 - (^^ip) « 1-32032 .... (6-12) 
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Physically the Hardy-Littlewood conjecture gives the number of prime pairs such 
that pi — P2 = r and pi < N 

N{p,p + r are primes and p < N) ^- — 5° ^ 2jy '^(^) " (6' 13) 

In particular it predicts that the number of twins primes (i.e. pi—p2 = 2) is asymp- 
totically C2N/ \p? N . This conjecture is well confirmed by existing numerics but no 
rigorous proof exists. The heuristic derivation of the Hardy-Littlewood conjecture is 
given e.g. in Ref. 14). 

§7. Two-point correlation function of Riemann zeros 

Taking the above formulae as granted one gets 



{off) 

■2 47^2 ^ „ 



After substitution the Hardy-Littlewood formula for a{r) (j6-8p and performing the 
sum over all r one obtains 

n (p,g)=l ^^^^^^ ^ 

Changing the sum over n to the integral gives 

2 



(p,g)=l ^^^'^^^ ^ 

The sum over coprime integers can be performed by using the inclusion-exclusion 
principle 

00 

E /(p) = EE/(^^H^)' M 

(p,g)=l k=l S\q 

and taking into account that 271(1 = ln(i?/27r). Finally one obtains''^ 

d°ff\e) = ^|C(1 + ie)|2e2-*'^"^cp(°//)(e) + c.c. . (7-5) 
The function •'^^ (e) is given by a convergent product over all primes 

From the above formulae the two-point correlation function of Riemann zeros is the 
sum of three terms 

R^[e) = d'iE) + Rf'''\e) + B!i^^\e) (7-7) 
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where the mean density of the Riemann zeros, d{E), is given in (j3-3p . the diagonal, 
R^^°'^\e), and ofF-diagonal, R^^'^\e), contributions are defined in (|6-3p and (j7-5p 
respectively. 

The very optimistic error in these formulae is supposed to be of the order of 
1 / which means that for all practical reasons they can be considered as the exact 
ones. 

To compare them with numerics it is necessary first to perform the unfolding 
of the spectrum. It is done by measuring the distance between zeros in units of the 
mean zeros density and re-scaling the correlation function as follows 

Here e is considered clS cl finite value. 

The limit S — > oo implies d{E) — > oo and the above formulae tend to the 
universal GUE result. 

2 — (27rg)2 ^ + ' ^ ' 

which agrees with oscillating terms of the GUE prediction ()4-5p . In Ref. 25) it was 
demonstrated that in this limit higher order correlation functions of the Riemann 
zeros also tend to the GUE results. 

The importance of the above formulae comes first of all from the fact that they 
describe not only the universal (structureless) GUE limit but also non-universal 
approach to it. 

Taking into account the first corrections to the GUE result one gets*^) 

R2{e) = 1 - ^^4^ - ^sin2(^e) - -i^^e sin(27re) + 0(J-^) . (7-10) 
Here /3 and 6 are numerical constants 

P = ll + 271 + E « 1-57314 , S = Y^ « 2.3157 (7-11) 

where the summation is performed over all prime numbers and 7„ are the Stieljes 
constants defined by the limit 

m^oo \ k m + 1 I 

From Fig.[2]it seems that numerical calculations agree well with the GUE prediction. 
Nevertheless, if one subtract the GUE result (|4-5p from these data it becomes evident 
that in the difference there is a clear structure (see Figs. [3]and[3]). In Fig. [5] (left) 
we present the difference between numerical two-pooint correlation function and the 
derived formula for it (|7-7p . The result is structureless and its histogram is given in 
Fig. [5] (right). Solid line in this figure represents the Guassian fit to the histogram. 
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Fig. 4. The same as in Fig. |3] but in the interval 10 < e < 15 (left) and 15 < e < 20 (right). 



The width of this Gaussian is close to the value of the bin used by Odlyzko in the 
computation of the two-point function. This comparison is the best verification that 
the theoretical formula (j7-7p extremely well described monumental calculations of 
Odlyzko. The difference is purely statistical and shows no visible structure. 

§8. Nearest-neighbor distribution 

There exist heuristic methods which permit the calculation of higher order cor- 
relation functions'^^-* but the resulting expressions become more and more tedious. 
For example, the three-point correlation function of Riemann zeros has the form^^ 

^3(61,62,63) =i?f"'^ (61,62,63) +4°''^ (61,62,63) (8-1) 

where the diagonal part is 
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Fig. 5. Left: the difFerence between the two-point function computed by Odlyzko and the theoretical 
prediction (|7-7p . Right: The histogram of the previous differences. Solid line is the Gaussian 
fit. 



+ 



1 



+ 



1 



(pl-ie2i _ l)(pl-ie23 _ 1) (pl-ie32 _ l)(pl-ie3i _ 1) 

and the oscillating part is given by the following formula 



+ c.c. 



.(off), 



61,62,63) 



(27ri)3 



iel2^2 



(p-l)2 



_d_ 



+ 



■log 



C(l + i632) 



c(i+i63i: 

ei2 _ 2 



o«ei2 



1 



(•pl+ie23 
1612^2 



l)(p 



l+iei3 



+ 



(pl-iei3 _ l)(pl-ie22 _ 1; 

-|-permutations + c.c. . 



+ 



(1 _piei2) 

p - 2 + p'ei2 



V 



1631 



+ 



1 



V 



l-ie23 



(8-3) 



Here eij = — ej and it is assumed that the cyclic permutations of indices 1, 2, 3 is 
taken together with the complex conjugation of the answer. 

The direct verification of these formulae is not easy but they can be used to 
obtain corrections to the other important statistics, namely, the nearest-neighbor 
distribution for Riemann zeros (cf. Fig. [2] (right)). This can be done as follows. 

For standard ensembles of random matrices correlation functions have the de- 
terminantal form (j4-3p . In many cases from physical considerations it follows that 
first order corrections to standard random matrix ensembles correspond to a change 
of the kernel K{Ei,Ej) (j4-4p only. In Ref. 8) this statement has been checked using 
the two-point and the three-point correlation functions presented above and it was 
conjectured that it remains true for higher correlation functions of Riemann zeros as 
well. 

From (j7T0p it is easy to check*^^ that the modified kernel takes the form 

K{e) = Ko{e) + h{e) (8-4) 
where -fCo(^) is the universal kernel (j4-4p and ki{e) is the effective correction to the 
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universal result 

where parameters A'gff and a have the following values 

_ TTdjE) \n{E/2-K) _ 5 _ 5 

VW ' "~ 27Td{E)(3~ ^ pin{E/27r) ■ 

Here P and (5 are the same constants as in (|7-lip . 

Using these expressions it was shown in Ref. 8) that the nearest-neighbor spacing 
distribution of the Riemann zeros can be calculated as follows. 

• First it is necessary to find the expansion into powers of N of the nearest- 
neighbor distribution for a circular unitary ensemble of x random unitary 
matrices equipped with the Haar measure (called CUEat) 

p^^---\s) = pois) + ^ pr"H^) + 0{N-') , (8-7) 

• Then the nearest-neighbor distribution for the Riemann zeros equals to the 
universal random matrix result, po{s) plus the correction terms Sp{s) where 

^Pi') = Pr''(«^) + OiN^s') ■ (8-8) 

with and a as above. 
Fig. shows the comparison between the difference of the numerical result of 
Odlyzko and asymptotic formula Eq. (j8-8p for a billion zeros located on a window 
around E = 2.504 x 10^^ (as in Fig. [2] right). The effective matrix size is N^s = 7.7376 
and a = 1.0438. The agreement is quite good. For comparison, we have plotted as 
a dashed curve the theoretical formula (j8-8p without the rescaling of the s variable. 

Fig. mb) is a plot of the difference between Odlyzko's results and the prediction 
(j8-8p . There is still some structure visible which might be attributed to the 0{N^) 
correction. To test the convergence, we have made the same plot but now using one 
billion zeros located on a window around E = 1.307 x 10^^ which corresponds to 
A^eff = 11.2976 and a = 1.0300 (Fig. El^a)). Now the agreement is clearly improved. 
The difference between the prediction (j8-8p and the numerical results, plotted in 
Fig. mb) shows a structureless remain. 

§9. Moments of the Riemann zeta function 

The purpose of this Section is to discuss the behavior of mean moments of the 
Riemann zeta function 

rT 



M,{T) = ^l^ |C(l/2 + it)pMt (9-1) 



when T oo. 



Riemann zeta function and quantum chaos 



15 




Fig. 6. (a) Difference between tfie nearest neiglibor spacing distribution of the Riemann 
zeros and the asymptotic GUE distribution for a biUion zeros located in a window near 
E — 2.504 X 10^^ (dots), compared to the theoretical prediction Eq. (|8-8|) (full line). 
The dashed line does not include the scaling of s. (b) Difference between the numerical 
Riemann values (dots) and the full curve (theory) of part (a). 



(a) 



(b) 



>5 



s s 

Fig. 7. Same as in Fig. 2 but for a billion zeros located in a window near E — 1.307 x 10^^. 



Naive calculations of this quantity can be performed as follows. Using the rep- 
resentation of C{s) as the Euler product (j2-2p one formahy gets 

|C(l/2 + it)P^ = J|(l - ylpe-^^i°P)-\l - Ape'^^'^P)-^ (9-2) 
p 

where Ap = p^^/^. 

Logarithms of prime numbers, \np, with different primes are non-commensurated. 
Therefore, for any finite number of primes e^-^''^^' plays the role of random phase e"^ 
and 

lim - / F{e'^^''P\e''^^''P\...e'^^''P'')dt= (9-3) 
^^oo T Jq 



T^oo T 

f2TT f2n f2TT 



JQ 



i^^e ,e ,...,e J 27r 27r " " " 27r ' 
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In application to ()9-2p it is necessary only to compute the integral over factors with 
the same p which can easily be done by the using the binomial expansion 

m=Q ^ ' 







The result is 

^"(1 - A,e--^^^)-\l - M^n-'^ = E i^^^y. (9-5) 

m=0 ^ ■ \ I / 

Finally, it is tempting to conclude that mean moments of the Riemann zeta function 
are the product over all prime of factors (j9-5p 



TTTT^MxiT) = a(A)/(A) (9-8) 



Unfortunately, at the critical line (i.e. when Ap = p ^/^) the product diverges and 
the divergent part is Yipi^ + ^ /p)- To get a finite expression let us define 

p m=0 ^ ■ \ I / 

The factor (1 — 1/p)'^^ ensures the convergence of the infinite product. 

Due to the divergence of the product in (j9-6p it was evident that certain increas- 
ing factors of T have to be introduced. So it was conjectured that instead of (j9-6p 
one should have the following result 

1 

tZ"oo (In'r)^ 

with an unknown function /(A). 

This conjecture has been based on the two proved theorems: /(O) = 1 by defini- 
tion, /(I) = 1 (Hardy and Littlewood, 1918), /(2) = 1/12 (Ingram, 1926), and two 
conjectures: /(3) = 42/9! (Conrey, Ghosh, 1992) and /(4) = 24024/16! (Conrey and 
Gonek, 1998). 

The breakthrough in this problem was the work of Keating and Snaith.^^ Their 
argumentation was as follows. By definition, C(l/2 + i-E'j) = where Ej are non- 
trivial zeros . Therefore, one can write C(l/2 + i-E) ~ nj(-^ ~ ^j)- such zeros are 
in a certain sense eigenvalues of an unitary matrix U this representation means that 
the zeta function is proportional to a characteristic polynomial of such matrix i.e. 
C(l/2 -|- \E) ~ det(e^'^ — U). But for random matrices from standard ensembles the 
moments of characteristic polynomials can be computed explicitly. 

In partricular for N x N unitary matrices equipped with the Haar measure the 
calculations are straighforward. By definition one gets 

2A 



U{N) 



^(nf^) n 

\k=l-^^ / l<j<m<N 



N 



n=l 



gi9 _ 



(9-9) 
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The integrals can be computed from the Selberg integral (see e.g. Ref. 21)) and the 
result is the following 



N 



2A\ _i-r r(i)r(2A + j) 



When N ^ oo the right-hand side of this expression tends to 



^P^)^^^/-^°°ArAVcoMA) (9-11) 

where the function fcOE{^) is 

G\l + \) 

= G(1 + 2A) ^^"^^^ 

and G(^;) is the Barnes G-function 

oo 

G(l + z) = (27r)^/2e-[(i+T)^'+^l/2 J] J(l + z/n)"e-^+^'/(2")] . (9-13) 

n=l 

By analogy with r{z) function the function G{z) can be defined as follows G(l) = 1 
and G{z + 1) = r{z)G{z). (Cf. r(l) = 1 and r(z + 1) = zr{z).) 

Let us compute the first values of fcOE{^)- fcOEi'O) = 1, /co£;(l) = 1) 
fcoEil) = 1/12, /coe(3) = 42/9!, /coe(4) = 24024/16!. Comparing with proved 
and conjectured values of /(A) cited above Keating and Snaith had conjectured that 
for general A 

/(A) = fcoE{\) (9-14) 
and the matrix dimension, A'^, asymptotically is related with T as 

Ar = lnr. (9-15) 

Under this identification all known and conjectured values of zeta function moments 
are in agreement with values obtained from the random matrix theory. 



§10. Conclusion 



• Quantum chaos methods and ideas are extremely successively applied to the 
Riemann zeta function and other zeta functions and L-functions from number 
theory. 

• Obtained heuristic formulae are very well confirmed by existing numerics. Even 
tiny details of monumental numerical calculations of Odlyzko are explained by 
these formulae. 

• There exists a quite large number of random matrix type conjectures for 'all' 
possible quantities like spectral statistics, moments, ratios of different zeta func- 
tions etc. 

• These conjectures are well accepted but no one of such results is proved rigor- 
ously. 
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